Dynamics of Granular Stratification 
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Spontaneous stratification in granular mixtures — i.e. the formation of alternating layers of small- 
rounded and large-faceted grains when one pours a random mixture of the two types of grains 
into a quasi-two dimensional vertical Hele-Shaw cell — has been recently reported by H. A. Makse 
et al. [Nature 386, 379 (1997)]. Here we study experimentally the dynamical processes leading 
to spontaneous stratification. We divide the process in three stages: (a) avalanche of grains and 
segregation in the rolling phase, (b) formation of the "kink" — an uphill wave at which grains are 
stopped — at the bottom substrate, and (c) uphill motion of the kink and formation of a pair of 
layers. Using a high-speed video camera, we study a rapid flow regime where the rolling grains 
size segregate during the avalanche due to the fact that small grains move downward in the rolling 
phase to form a sublayer of small rolling grains underneath a sublayer of large rolling grains. This 
dynamical segregation process — known as "kinematic sieving", "free surface segregation" or simple 
"percolation" — contributes to the spontaneous stratification of grains in the case of thick flows. We 
characterize the dynamical process of stratification by measuring all relevant quantities: the velocity 
of the rolling grains, the velocity of the kink, and the wavelength of the layers. We also measure 
other phenomenological constants such as the rate of collision between rolling and static grains, and 
all the angles of repose characterizing the mixture. The wavelength of the layers behaves linearly 
with the thickness of the layer of rolling grains (i.e., with the flow rate), in agreement with theoretical 
predictions. The velocity profile of the grains in the rolling phase is a linear function of the position 
of the grains along the moving layer, which implies a linear relation between the mean velocity and 
the thickness of the rolling phase. We also find that the speed of the upward-moving kink has the 
same value as the mean speed of the downward- moving grains. We measure the shape and size of 
the kink, as well as the profiles of the rolling and static phases of grains, and find agreement with 
recent theoretical predictions. 
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I. INTRODUCTION 

Size segregation of granular mixtures [j]^]] is known 
to occur when mixtures are exposed to external peri- 
odic perturbations. A much-studied size segregation phe- 
nomenon is known as the "Brazil nut effect" [jsl-fi^f and 
occurs when, upon vibration, larger grains rise on a bed 
of finer grains. Axial size segregation in alternating bands 
consisting of small and large grains occurs when a mix- 
ture of grains is placed in a horizontal rotating cylinder 
p3| |r6[ . ft is also known that even in the absence of ex- 
ternal perturbations mixtures of grains of different sizes 
can spontaneously segregate. For example, when a mix- 
ture of spherical grains of different sizes is poured onto 
a heap, the large grains are more likely to be found near 
the base, while the small grains are more likely to be near 
the top 

Another type of segregation, called spontaneous strat- 
ification, arises when the grains composing the mixture 
differ not only in size but also in shape (or friction prop- 
erties). When a mixture of large grains that are more 
faceted and small grains that are less faceted is poured 
in a "granular Hele-Shaw cell" (two vertical slabs sepa- 
rated by a gap of typically 5-10 mm), the mixture sponta- 
neously stratifies into alternating layers of larger faceted 
grains and smaller rounded grains [^7j . Figure [l]a shows 
an example of such stratification. A mixture of large cu- 



bic sugar grains (typical diameter 0.8 mm) and smaller 
spherical glass beads (diameter 0.19 mm) is poured in the 
cell. We notice the striped pattern with approximately 
constant wavelength. 

In contrast, when the mixture is composed of larger less 
faceted grains and smaller more faceted grains, the mix- 
ture only segregates — i.e., the small more- faceted grains 
are found preferentially at the top of the cell, while the 
large less-faceted grains are found near the bottom [p7| . 
Figure [l]b shows an example of such segregation, when 
a mixture of small faceted sand grains (typical size 0.3 
mm) and large spherical glass beads (typical size 0.8 mm) 
is poured in the cell. 

The dynamical process leading to spontaneous strat- 
ification was recently studied theoretically [28|-|3l|l us- 
ing discrete models, and a set of continuum equations 
for surfaces flows of granular mixtures developed in 
Refs. [^2 3 5 1 . In this theoretical formalism, the grains 
are considered to belong to one of two phases: a static or 
bulk phase if the grain is part of the solid sandpile, and 
a rolling or liquid phase if the grain is not part of the 
sandpile but rolls downward on top of the static phase. 
In Ref. |28| the dynamics of spontaneous stratification 
was found to be governed by the existence of a "kink" at 
which the grains are stopped during an avalanche. 

In this paper, we study experimentally the dynamical 
processes leading to spontaneous stratification. Using 
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a high speed video camera to study the motion of the 
grains in great detail, we divide the dynamical process of 
stratification into three stages (see Fig. ||): 

(a) The avalanche of grains down the slope, and size 
segregation of grains in the rolling phase due to 
"percolation" . 

(b) The formation of the "kink" — an uphill wave at 
which grains are stopped. 

(c) The uphill motion of the kink and formation of a 
pair of layers. 





(b) 

FIG. 1. (a) Example of stratification for a mixture of 
smaller rounded grains (white spherical glass beads) and 
larger faceted grains (black sugar grains), (b) Example of 
sharp segregation for a mixture of smaller faceted grains 
(black sand) and larger rounded grains (white spherical glass 
beads). Notice the sharp zone of separation of the species in 
the middle of the pile. This is the result of strong segregation 
effects acting in the system. Notice also the smaller angle of 
repose of the spherical beads at the bottom of the pile. 

We study a well-developed flow regime where the 
rolling grains segregate during the flow. In this regime 
the thickness of the layer of rolling grains is larger than 



the typical size of a grain d (typically 5d) , and the smaller 
rolling grains are found to percolate downward in the 
rolling phase to form a sublayer of smaller rolling grains 
underneath the sublayer of larger rolling grains. This dy- 
namical size segregation process , known as "percolation" 

, contributes to the 
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or "kinematic sieving" 
stratification of grains. 

Stratification is an instability developed due to a com- 
petition between size segregation and shape segregation 
p9[ . In the case of thin flows, size segregation occurs 
since the smaller grains are captured more easily than 
larger grains. In the case of thick flow regimes study here, 
the kinematic sieving in the rolling phase is mainly re- 
sponsible for the size segregation of the grains. Since the 
larger grains are on top of the rolling phase, they are con- 
verted further down than the smaller grains, producing 
the size segregation effect, which together with the seg- 
regation due to different shape of the grains, gives rise to 
the instability leading the system to spontaneously strat- 
ify p9[ . It is important to note that percolation in the 
rolling phase is not sufficient condition to obtain strati- 
fication. For thick flows and when the large grains are 
smoother, segregation in the rolling phase still occurs, 
and yet we do not get stratification but only the sharp 
segregation pattern of Fig. [l]b. 

A large difference in size is also a condition for the per- 
colation effect to take place — usually p > 1.5, where p 
is the ratio of the size of the large grains to the size of 
the small grains. We performed a series of experiments 
with mixtures of glass beads and sand with p < 1.5 and 
found continuos segregation patterns (as opposed to the 
sharp segregation pattern with a separation zone of a few 
centimeters of Fig. [I|d obtained for p > 1.5) no matter 
the shape of the grains. This is because, when p < 1.5 
the effect of size segregation is very weak. 

The limiting case in which both species of grains are 
spherical was first studied by Williams [^9| , ^6|j37| ; his re- 
sults (showing segregation plus a hint of stratification) 
differ from our results for this case (showing only seg- 
regation). We believe that the reason is that the grains 
used by Williams were not quite spherical, inducing some 
shape segregation as well. According to the above in- 
terpretation, we note that experiments with mixtures of 
perfect spherical beads differing only in size should not 
show stratification since the shape segregation effect is 
not present and size segregation alone (even due to per- 
colation) is not able to produce stratification — and our 
work confirms these expectations. However some oscilla- 
tions might still be present around the stable segregation 
profile, as seen in previous experiments using mixtures of 
spherical beads Q . 

Here we focus on the regime where segregation in the 
rolling phase takes place. We characterize the kinematic 
percolation process, and measure the velocity gradient 
of the grains during the avalanche. We find a linear ve- 
locity profile of the rolling grains, and that the mean 
velocity of the rolling grains is the same as the velocity 
of the kink. We also study the shape and size of the 
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kink — and thereby measure the wavelength of the layers. 
We find that the wavelength increases linearly with the 
flux of grains, a result in agreement with recent theoret- 
ical predictions |28| . We also measure the profiles of the 
rolling grains and static phases of the pile, and the values 
of several phenomenological coefficients which appear in 
the theory for surface flows of granular mixtures. Our re- 
sults are valid for flow rates of the order of gr/sec (which 
gives rise to a rolling phase less than 1 cm thick). We 
also comment on the applicability of our results, and on 
the deviations that may occur for smaller and larger flow 
regimes. 




FIG. 2. Three stages of the dynamics of stratification: (a) 
avalanche of grains, (b) formation of the kink at the bottom 
substrate, (c) uphill motion of the kink. The dashed line in 
(b) is the window of observation used to record the images. 



II. EXPERIMENTAL SETUP 

Our experimental setup consists of a granular Hele- 
Shaw cell: a vertical "quasi-two-dimensional" cell with a 
narrow gap separating two transparent plates (made of 
plexiglass, or of glass). The cell measures L = 30 cm of 
lateral size and 20 cm high, and the gap is I = 0.5 cm. 
We close the left edge of the cell. We clean the walls of 
the cell with an antistatic cleaner in order to avoid the 
effects of electrostatic interaction between the grains and 
the wall. 

In this study, we focus on spontaneous stratification. 
In all our experiments, we use a mixture of grains com- 
posed of two species differing in size and shape: smaller 
glass beads of average diameter 0.19±0.05 mm, spherical 
shape (95% sphericity), angle of repose On = 26° ± 1° 
(we call these type 1 grains), and larger faceted sugar 
grains of typical size 0.8 mm, approximate cubic shape 
and angle of repose 622 = 39° ± 1° (type 2 grains). 

The typical size of the sugar grains obtained by mea- 
suring the volume of the cubic grains and calculating the 
typical size as the cubic root averaged over 20 different 
grains. We obtain the value of the angle of repose of the 
species by pouring the pure species in the Hele-Shaw cell 
and measuring the resulting angle of the pile, averaging 
over 5 realizations of the sandpile. The angle thus mea- 
sured is not the actual angle of repose corresponding to 
a conical pile, since the presence of the wall induces ex- 
tra friction that slightly increases the equilibrium angle 
of the pile However, we are interested in the angle 
of repose for this specific geometry since our experiments 
on stratification are done in the cell. 

We fill the cell at different rates of adding grains (flux) . 
We use a Kodak Ektapro 1000 digital high speed camera 
to film the motion of individual grains during the for- 
mation of the layers. The camera produces 1000 digital 
frames per second with a resolution of 238 x 191 pix- 
els. We record images during 1.6 sec, and achieve longer 
recording times by lowering the frame rate. The digital 
images we download to a workstation for further image 
processing. 



III. THE ANGLE OF REPOSE OF THE PURE 
SPECIES 

Since stratification is related to the different angles of 
repose of the pure species, we first study how the angle 
of repose depends on the size and shape of the grains. 
We measure the angle of repose of different sets of spher- 
ical glass beads of different size, and find the same value 
within errors bars (see Fig. ||). The angle of repose 
does not depend on the size of the grains, since a simple 
isotropic rescaling of the pile coordinates transforms a 
pile of smaller spherical grains into a pile of large spher- 
ical grains, while leaving invariant the angle of the pile. 
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FIG. 3. Angle of repose of different set of spherical glass 
beads differing only in size, and other species such as sand, 
salt and sugar. The more irregular is the shape of the grains, 
the larger is the angle of repose. 

The value of the angle of repose we find for the spher- 
ical glass beads is smaller than the value we find for the 
cubic sugar grains. In general, we find that the angle of 
repose does not depend on the size of the grains, and is a 
function of the shape of the grains: the rougher the shape 
of the grains the larger the angle of repose — because for 
more faceted grains the packing of grains is less dense 
than for more rounded grains. 

A particular case is found when the size of the grains 
is of the order of microns. Spherical beads of size 40/im 
have a larger angle of repose than the same spherical 
beads but of size of the order of mm (see Fig. ||) — because 
adhesion forces become important, increasing the angle 
of repose of the species. The scale of microns is the lower 
limit of applicability of our results, since at the submi- 
cron scale particles undergo Brownian motion [Q and our 
analysis of collisions and transport at zero temperature 
ceases to be valid. 



IV. THE KINK MECHANISM 

A physical mechanism has been proposed for the for- 
mation of the layers which is related to the existence of 
a "kink" p7p^ , |30| ]. Suppose, e.g., that a pair of static 
layers is formed with the layer of larger grains on top of 
the layer of smaller grains. When an incoming mass of 
grains avalanches down the slope, the larger grains reach 
the base of the pile first, due to the fact that large grains 
do not tend to get trapped (in local minima of the sand- 
pile profile) as easily as small grains. Additionally, in 
the case of rapid flows, the smaller grains also size segre- 
gate to the bottom of the rolling phase due to percolation 
so that the larger grains, being at the top of the rolling 
phase, tend to travel further since they do not interact 
with the bulk phase. 

During the avalanche of grains, some small grains are 
captured in the static layer of large grains, smoothing 
the surface and thereby allowing more small grains to 



fall downward and eventually reach the bottom of the 
pile. When the flow reaches the base of the pile, we 
see that the grains develop a profile characterized by a 
well-defined "kink" at which the grains are stopped. This 
kink moves in the direction opposite to the flow of grains, 
conserving its profile until it reaches the top of the pile. 

In the process of falling down the slope, grains (small 
and large) stop at the kink. We see that the smaller 
grains stop first (since the small grains are already segre- 
gated in the rolling phase) so a pair of layers forms, with 
the smaller grains underneath the large grains (see Fig. 
|l|b) . When the kink profile reaches the top of the sand- 
pile, the pair of layers is completed. Then this process 
is repeated: a new avalanche occurs, some larger grains 
reach the bottom of the pile, the kink is developed, and 
a new pair of layers is formed. 

The size of a pair of layers A is determined by the thick- 
ness of the layer of rolling grains during an avalanche, R , 
which in turn is determined by the flux of adding grains. 
The volume of rolling grains O ava i that reaches the kink 
during a time interval At and in a differential dy is 



aval 



£At v(y)dy + £At v^dy, 



(4.1) 



where v(y) is the velocity of the rolling grains at a dis- 
tance y from the pile surface of static grains, > is 
the upward velocity of the kink which is constant, and £ 
i s th e gap between the plates of the cell. The first term in 
(4.1) is the volume of grains falling down the slope, and 
the second term represents the volume of grains from the 
rolling phase that the kink encounters when it advances 
uphill at velocity v\ — i.e., v(y) +i>i is the velocity of the 
rolling grains in the co-moving frame of the kink. 
The volume of grains in a well-formed kink is 



f^kink = £At v^dy. 



(4.2) 



Then if all the grains are stopped at the kink, the num- 
ber of rolling grains falling down, ^tfiuid^avai, where /Zfl u id 
is the density of the fluid phase (the number of rolling 
grains per unit volume) should scale as the volume of 
grains in the kink /Ubulk^kinkj where /Xbuik is the density 
of the bulk phase. Hence 



(v(y) + vf) dy = fi huik £At / v^dy. 



(4.3) 



The mean value of the downward velocity of the grains 
averaged over the rolling phase is 



i r 

v = -j$; J a v (v) d y> 

so from (O) |7]J2| 



A = ^fluid {V + Vf) r0 



Mbulk 



,; 1 



(4.4) 



(4.5) 
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The analytical shape of the kink has been obtained in 
|{28| , j30f . We introduce four different generalized angles of 
repose a p, corresponding to the interactions between a 
rolling grain of type a and a static grain of type (3: 

• 6*22 corresponds to the angle of repose of the pure 
large-cubic species, 

• 6*ii is the repose angle of the pure small-rounded 
species (6* 2 2 > and 

• 621 correspond to the interaction between a large 
cubic rolling grain and small-rounded static grains, 
and 

• #12 correspond to the interaction between a small- 
rounded rolling grain and large cubic static grains. 



For stratification we have u 

O21 < 0\\ < $22 < $12 



(4.6) 



Since the kink is a traveling wave solution, we can write 

IS 



grains and therefore has an angle close to On . At the cen- 
ter of the kink, the larger grains are captured on top of 
smaller grains, therefore the angle decreases toward the 
cross-angle of repose #21 ■ Then the angle of the profile of 
the kink approaches 8q. 

Dimensional analysis of Eq. (T8) indicates that the 
upward velocity of the kink is proportional to the flux of 
arriving grains, i.e. proportional to R° [E8| , 



Ci 7 R°- 



(4.12) 



Here C\ is a numerical constant that does not depend 
on 7 or R°, but may depend on the angles of repose and 
other features of the grains. 

Below, we will show that the velocity of the rolling 
grains in an avalanche increases linearly with the height 
of the grains in the rolling phase, implying that the mean 
value of the velocity v is proportional to the thickness of 
the rolling phase. The coefficient is again proportional 
to 7 by dimensional analysis, 



f(x, t) = h(x, t) + 6 Q x = f(u), 



(4.7) 



C 2 7 



(4.13) 



where u = x + v ft, 9q is the angle of the pile after a pair 
of layers is formed, and h(x, t) is the height of the static 
phase (Fi g. [2p ) . The solution for the lower layer of small 
grains is pq] 



1 



log 1 



2wf 
R° 



1_ 

V] 



(/ - S lU ), 



where 



61 = O - On > 0, 



(4.8) 



(4.9) 



w = Uf /(v + Uf ), and 7 (units of 1/ sec) is the rate of col- 
lisions between static and rolling grains. Since we have a 
pile made of two different species, the angle of the layers 
is not necessarily either 622 nor On. However, since 
the top layer of the stripes is made of large-cubic grains, 
the resulting angle 0q is closer to #22 than to On. 

The solution for the profile of the upper layer of the 
kink is il 



R° 

f(u) = [ - 
w 



1 — exp 



where 



S2 = 0o - 1 



'22 



< 0. 



(4.10) 



(4.11) 



We will compare this solution with the profile of the kink 
measured experimentally. Figure |^c shows a sketch of 
the kink. The angle of the layers is #0, which is an angle 
between 9\\ and O22, so that solutions ( [1.8[ ) and (4.10) ex- 
ist. The lower part of the kink is made of small-rounded 



where C2 is also a numerical constant that does not de- 
pend on 7 or R , but may depend on the angles of repose. 

From Eqs. 



5|), ([yj) and 13| ) we obtain the de- 
pendence of the wavelength on R® 



Mfiu 



Mbulk 



£ CR°, 



(4.14) 



where C = 1 + C^/Cx, is a constant independent of R°. 
Thus, the wavelength incr ease s linea rly with the flux of 
grains. We will test Eqs. (4.8)-( 4.14 ) experimentally. 
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FIG. 4. Trajectories (x,y) of the rolling grains during the 
avalanche for Exp. #5 in a window of observation of 2.82 cm 
by 2.26 cm. The thickness of the rolling layer is R° = 0.65 
cm. The dashed line is the boundary between the sublayer of 
smaller rolling grains and the sublayer of larger rolling grains 
due to the percolation effect. The angle of the pile is 60 ■ 
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V. DYNAMICS OF STRATIFICATION 

To test the above mechanism, we perform a series of 
six experiments where we study in detail the dynamics of 
stratification by measuring all the quantities involved in 
the process. We use a Kodak Ektapro 1000 high-speed 
video camera system to record the motion of the grains. 
In order to study the profile of the kink and the effects 
of percolation in the rolling phase, we measure the pro- 
file of the static and rolling phases and compare with 
analytical predictions. We measure the velocity of the 
rolling grains, the velocity of the kink, the wavelength 
of the layers, and also other phcnomcnological constants 
such as 7, and all four angles of repose characterizing the 
mixture 9 a p. 



According to the picture discussed in Sec. IV, we study 
the dynamical process of stratification by dividing the 
process in three stages (Fig. ||): 




FIG. 5. (a) Position x(t) and (b) position y(t) of the rolling 
grains shown in Fig. ^ as a function of time. The straight lines 
indicate that the velocities achieve a constant value, given by 
the viscous damping and the gravitational driving force. The 
symbols correspond to the grains plotted in Fig. 



We focus our study on a small window of observation 
of typically 3 cm x 2 cm size located in the center of 
the pile. Using the high speed video camera at a frame 
rate of 1000 frames per second, we are able to track the 
motion of each individual grain in a downhill avalanche 
during 1.6s of recording time. The rolling grains take 
only 0.1s to cross the window of observation. However, 
we continue recording after this time elapses, so that we 
can record in the same shot of 1.6s the grains flowing 
down and the profile of the kink moving up. By tracking 
the motion of each individual grain, we are able to mea- 
sure the velocity profile of the grains along the layer of 
moving grains. 

The thickness of the layer of rolling grains in all our 
experiments ranges from 0.3 cm to 0.7 cm. Thus the 
layer of rolling grains is thick enough that it is possible 
to observe the size segregation of the grains in the rolling 
phase. For a well-developed flow of grains down an in- 
clined plane, the grains segregate in the moving layer, 
with the small grains at the bottom of the moving layer, 
and the large grains at the top. This effect is called "free 
surface segregation", "percolation" or "kinematic siev- 
ing" |^,|l],|3]-|2^] , and occurs because the smaller grains 
percolate downward through the gaps left by the motion 
of larger grains in the rolling phase. In our experiments, 
the flux of grains is sufficiently high that the layer of 
rolling grains is large enough that the percolation effect 
can be observed in the rolling phase. Figure ^ shows the 
trajectories of several grains for Experiment #5 (Table 
||) in a window of 2.82 cm by 2.26 cm, plotted every two 
frames (« 2ms). We find that the large grains occupy 
the top part of the moving layer, and the small grains 
are located at the bottom of the moving layer. 
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FIG. 6. Velocity parallel to the pile surface of the rolling 
grains, corresponding to the experiments listed in Table |. as 
a function of their vertical position y from the top of the static 
phase. We find a linear velocity profile. 



A. Avalanche of grains and percolation effect 

In all six experiments, we use the same mixture and 
cell, but different fluxes — i.e., different BP (see Table |). 



To measure the velocity profile of the grains in the 
rolling phase, the position of every rolling grain is spot- 
ted on the screen of the video camera. We follow the 
trajectory of the grain during a period of time where the 
rolling grain is well distinguished from the other grains 
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in the moving phase. We stop the image every 2ms (2 
frames) and record the (x, y) position in pixels of the 
screen. The position thus measured is manually entered 
in a data file, giving the (x, y) position of the grain as a 
function of time. We study the motion of the grains in 
the center of the pile, where the grains have achieved a 
constant velocity along the direction of the pile surface 
(the viscous friction force has balanced the gravitational 
driving force on the grain) . Figure ^ shows the x and y 
coordinates of the grains of Fig. |as a function of time. 
The velocity of the grains we calculate from the slope 
of such curves. Indeed, we observe that the motion of 
the grains in the center of the pile is overdamped; the 
velocity is constant as a function of time. 

Due to the percolation effect, the layer of rolling grains 
is actually made of two equal size sublayers (since we use 
an equal volume mixture of two species) of smaller and 
larger grains. However, the velocity profile of the grains 
is continuous along the thickness of the moving layer. 
Figure ^ shows the velocity of the grains for the six ex- 
periments listed in Table | as a function of y, the height 
of the rolling grain to the top of the static phase. The 
velocity profile is linear in y. Using the data from the 
experiments listed in Table |l] we find 



v ( y ) = (46/sec±2) y. 



(5.1) 



The coefficient of the linear relation (5.1) is independent 
of imp lyin g that the mean value of the velocity of 
the grains (4.4) is proportional to R° 



(23/sec±2) R°. 



(5.2) 



By comparing with Eq. ( 4.13| ) and assuming that the co- 
efficient C\ is of the order of one, we obtain an estimate 
of the rate 7 ~ 23/sec. Similar velocity profiles have 
been found in Ref. j40[ j , although these results (which 
were obtained for single-species grains falling down in- 
clined planes at different angles above the angle of re- 
pose) do not apply to our system, since we are interested 
in the velocity profile of the grains avalanching on a sur- 
face composed of large grains at the angle of the layers, 
as occurs in the stratification experiment. 



B. Formation of the kink 



The formation of the kink is determined by the inter- 
action of the grains at the bottom of the pile. When the 
grains reach the substrate, we find that the larger grains 
arrive first and then the smaller grains arrive (Fig. ^i), 
because the larger grains roll down more easily than the 
smaller grains, since they are at the top of the rolling 
phase due to percolation. They have larger velocity than 
the small grains since they are at the top of the rolling 
layer (Fig. |^a). The large grains are stopped at the sub- 
strate, and they develop the profile shown in Fig. (?]b. 

The condition for the formation of the kink seems to 
be the existence of two species, not one. In fact, when we 
pour single-species grains in a open cell, we do not ob- 
serve the stationary kink, but we observe that the height 
of the pile and the profile of the rolling phase acquire a 
steady state (without oscillations) where the profiles are 
conserved in time. In the case of two species, the larger 
grains reaching the bottom before the smaller grains act 
as a "wall" or "incipient kink" where the smaller grains 
are stopped (Fig. |^b). Thus, when the small grains ar- 
rive near the substrate, they find some large grains al- 
ready there. They are stopped in this way, and the kink 
emerges from this interaction. When the kink is devel- 
oped it starts to move uphill with constant velocity and 
conserving its profile. 

de Gennes has shown that when a flow of single- 
species grains flowing down a plane at the angle of repose 
reaches a vertical wall, the grains develop an uphill wave 
of constant velocity Uj ~ jR°. Although this uphill wave 
is not stationary as found for the kink in our experiment, 
the solution found by de Gennes shows that it is possible 
for the smaller grains to be stopped by a moving "wall" 
of large grains and thereby give rise to a kink. The ex- 
istence of uphill waves (although not stationary as the 
kink solution) was also noticed in Refs. ^J^f . 



C. Uphill motion of the kink and formation of a pair 
of layers 




FIG. 7. Formation of the kink at the base of the pile, (a) 
The large grains roll down ahead of the smaller grains since 
they are at the top of the rolling layer and have larger veloc- 
ity than the small grains, (b) When the larger grains reach 
the substrate they are stopped, and they act as a wall where 
the smaller grains are stopped. The kink emerges from this 
interaction. 



By using the video camera at a frame rate of 1000 
frames/sec, we can distinguish the fraction of grains 
which is in the rolling phase and the fraction of grains 
which is strictly immobile, the fundamental ingredient of 
the theories of Refs. §|-[||[^|3(|§|] . Thus we identify 
the time behavior of the boundary between the rolling 
phase and the static phase, and the profile of the kink. 

Since the contrast between the rolling and static grains 
is not very good, we must identify this boundary "manu- 
ally." We play the movie 5 frames back and then 5 frames 
forward and identify which grains are moving and which 
grains are static. We repeat these measurements every 
0.05s, and record the coordinates of the bulk/fluid inter- 
face 01. 
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Figure ||a shows the profile of the kink as a function of 
x plotted every 0.05 sec in a window of 2.82 cm by 2.26 
cm (Experiment #5), and one sees the profile of the kink 
moving upward with constant velocity v\. Notice that 
the profile of the kink is stationary. Figure |^b shows the 
angle of the kink profile 



9(x,t) = —axcta,n(dh(x,t)/dx), 



(5.3) 



where h(x, t) is the profile of the kink obtained from Fig. 
pk, as a function of x measured at different times. We 



fit th e analytical solution of the shape of the kink ( 4 
(4.10) obtained in pl| to the experimental profile of the 
kink and find good agreement (Fig. ||a). We find the best 
fit for values 7 = 150/sec, v = 15.6 cm/sec, v-\ = 16.2 
cm/sec, R° = 0.7 cm, 8 = 33°, 9 n = 25°, and 6» 22 = 43°. 




(a) 



1.00 2.00 
x (cm) 




(b) 



FIG. 8. (a) Profile of the kink h(x, t) obtained in Experi- 
ment #5, shown at time intervals of 50ms. The dotted line is 
the boundary between the static layers of smaller and larger 
grains. The kink moves uphill with constant velocity Vf . The 



symb ols co rrespond to the fits to the analytical solution (4.5) 
and ( |4.10| ). The data are first obtained manually and the 
fitted with a polynomial, (b) Angular profile of the kink ob- 
tained from Fig. (a), 6(x,i) — — arctan(9/i(i, t)Jdx). The 
angle is between 80 and Q21 as explained in Fig. Be. 



The solution (4.8) and (4.10) is valid only for the lower 
and upper part of the kink, so the center of the kink, 
where small and large grains are mixed, cannot be repro- 



duced. The values we use to fit the analytical solution 
are somehow different from the values we measure. How- 
ever, we note that the exact shape of the kink depends 
on the type of interaction term used in the equation of 
motion to describe the rolling/static grains interaction. 
In particular, the interaction term used in [p8[ does not 
include nonlinear terms that might be important when 
the flux of grains becomes large (we comment on this 
point in Sec. [v|). However, the fair agreement between 
experiment and theory indicates that some features of 
the interaction have been captured by this approach. 

We also focus on the different collision processes be- 
tween rolling grains in contact with the solid surface and 
the static grains. We find that amplification process (i.e. 
when a rolling grain removes a static grain via a colli- 
sion) do not happen very often. The percolation effect 
forbids the larger grains to interact with the surface, thus 
prohibiting cross-amplification processes of the type of 
a larger rolling grain amplifying a smaller static grain. 
The main interaction seems to be the capture of rolling 
grains at the surface — when a rolling grain is converted 
to the static phase. However, we emphasize that it is dif- 
ficult to clearly determine the nature of the interaction 
at the surface (capture versus amplification) because the 
smaller grains are the only interacting grains in the re- 
gion of observation and they are difficult to track. 

We measure the velocity of the kink Uf as a function of 
R° (see Table Q) . Figure ^a shows the results which can 
be fit to a straight line. We find 



j'T 



(23/sec±2) R°, 



(5.4) 



where R° is typically 5 cm. The velocity of the kink is ap- 
proximately the same as the mean veloc ity o f the rolling 
grains (C2 — C\). Comparing with Eq. (4.12), we obtain 
a second estimate of the rate 7 ~ 23/sec. 

Figure shows the wavelength of the layers as a func- 
tion of the thickness of the rolling phase for the six ex- 
periments of Table |. The data can be fit to a straight 
line, and we find 



A = (1.7 ±0.1) R°, 



(5.5) 



which agrees with the prediction of Eq. ( |4.14j) . Compar- 
ing with Eq. ( 4.14 ), we obtain (/Zfl u id/ Mbuik)C — 1-7. Us- 
ing (5.2) and (5.4) we obtain C ~ 2, so that /Xfl u id//J.buik — 
0.85, corresponding to the fact that the fluid phase is less 
dense than the bulk. 



VI. PHENOMENOLOGY 

Table |0| shows the values of the phenomenological con- 
stants measured for the equal volume mixture of quasi- 
spherical glass beads of mean diameter 0.19 mm and 
cubic-shaped sugar grains of typical size 0.8 mm. The 
mean value of the velocity of the grains falling down the 
slope and the velocity of the kink range from 7 cm/sec 
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to 17 cm/sec. As noted above, 7 « 23/sec and R° « 0.5 
cm. 



and 



57°. 



(6.3) 



0J 



< 

> 




(a) 



0.2 0.4 
R (cm) 



5 




R° {U) 



(b) 

FIG. 9. (a) Velocity of the kink as a function of the 
thickness of the rolling layer _R° for the six experiments (Ta- 
ble [i]). (b) Wavelength of the layers as a function of the 
thickness of the rolling layer R° for the experiments (Table 



According to Eq. ( |4.8j ) , the lower part of the kink has 
the angle of the small-rounded grains 9n, and near the 
layer of large grains there are logarithmic corrections. 
The angle of the kink decreases towards the center of the 
kink; it is equal to 621 at the position where the larger 
grains start to be captured (Fig. ^|c) Thus, measuring the 
minimum of the angle of the profile of the kink at the 
transition from the layer of smaller grains to the layer of 
larger grains provides a method to measure the value of 
the crossover angle of repose 6*21- Then, assuming that 
#ii — #21 = #22 — #12 p| , the critical angle 6x2 can be 
obtained too. 

We measured for the angle of repose of the pure species 

6> n = 26°±l, 6>22 = 39°±1. (6.1) 

Figure ||b shows the angle of the pile near the kink at 
different times. From these curves we measure the angle 
of the layers 

O ~ 33°, (6.2) 

and the remaining angle of repose of larger grains on 
smaller grains 



The constant u/7 represents the distance at which a 
rolling grain stops in a pile at an angle different from the 
angle of repose — i.e., v/j represents the distance at 
which a rolling grain is stopped at the kink. The con- 
stant v/j ~ 0.3 — 0.7 cm is small compared to the sys- 
tem size L = 30 cm — i.e., v/j scales with the size of the 
grains and not with L, as expected Q. Notice also that 
TJ/(7tan-0) — 2cm — where tamp = tan#n — tan #21 — 
0.3 — is the size of the region of mixing in the center of 
the pile observed in the case of segregation of the mix- 
tures of smaller cubic grains and large rounded grains 
(Fig. 0b) H^] . This region of the mixture is observed to 
be small in comparison with the systems size. 

Finally, we comment on the application of the theo- 
retical calculation of the model developed in p8| , p0| , p9]| 
to the problem of stratification when percolation effects 
are acting. The dependence of the re pos e angle on the 
composition of the surface proposed in p8| , ^o| , p9[ is analo- 
gous to the effect of percolation. Due to percolation, only 
the smaller grains interact with the surface, causing the 
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larger ones to be convcctcd further. The models of 
use the fact that the repose angle of the smaller grains is 
always larger than the repose angle of the larger grains 
for a given composition of the surface (i.e., 6\\ > #21, 
and 9i2 > #22)1 then the smaller grains are always the 
first to be trapped, and the large ones are always con- 
verted down as it occurs due to the percolation effect. 
Moreover, capture of larger grains is forbidden on top 
of smaller grains since the capture function of the large 
grains is zero around the angle of repose of the small 
grains |2S|j . 

A simple extension to ex plicitly include percolation ef- 
fects in the formalism of f2^,|3(| shows only small cor- 
rections to the profiles of the rolling and static grains, 
which provides evidence for the possible applicability of 
the results of 1 28 3^] to the case where percolation effects 
take place in the rolling phase. However, caution must 
be taken in the definition of the fluid/bulk interaction 
in the theoretical formalism. We take the interaction 
term to be proportional to the thickness of the rolling 
phase, an approximation valid for thin flows |7],|34|. Al- 
though this approximation might be still valid in the case 
of thick flows |Q — since the interaction might be pro- 
portional to the preassure exerted by the fluid phase [[13| , 
which in turn is proportional to R° for a fluid — nonlinear 
terms might be also necessary to completely describe the 
flow, especially because the interaction among the rolling 
grains (which is neglected in the theoretical formalism) 
becomes important. 

The dependence of the velocity of the rolling grains on 
the position of the grain in the rolling phase is another 
fact not included in the theory, which considers a uni- 
form velocity for all the grains in the rolling phase. For a 
comparison with the theory, we have replaced the veloc- 
ity of the grains in the theory by the mean value of the 
velocity of the rolling grains measured experimentally 
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VII. OTHER LIMITS 



VIII. DISCUSSION 



We notice that the dependence on the plate separation 
£, although pr esent in Eqs. (4.1) and (p~2]), has disap- 
peared in Eq. ( 4.14 ), and the relevant length scale that 
determines the wavelength is R . However, a change in 
the flux of adding grains, or a change in the gap i, changes 
the wavelength since R° ~ flux/^. Thus, e.g., by keeping 
the flux constant and increasing the gap I, we find a de- 
crease in R°, and w e find that the wavelength A decreases 
according to (4.14) 



A ~ l/l. 



(7.1) 



This dependence has been measured in J4^|43j . 

For the moderate fluxes used in this study (of the order 
of 1 gr/sec), the role of the flo w rate is to determine the 
wavel engt h according to Eq. ( 4.14 ). For larger fluxes, 
Eq. (4.5) is still valid as long as the kink mechanism 
works. Howe ver, deviations from the linear dependence 
of Eq. (4.14) might occur since the velocity of the kink 
and the velocity profile of the rolling grains might de- 
viate from the linear regime. The densities of the bulk 
and fluid phase might also change with the flux of grains, 
giving rise to nonlinear relation between A and R . For 
sufficiently large flow rates, the kink mechanism required 
to form layers cannot occur (especially the appearance of 
the first kink at the onset of the instability |2Sj] ), since 
the grains acquire large momentum, and avalanches that 
terminate by an upward moving kink before the next 
avalanche begins cannot occur. In this case, the kink 
is not able to stop the arriving rolling grains anymore; 
the grains ride over the kink so that no segregation at 
the kink is possible. Therefore, for this limit, the strat- 
ification pattern disappears when the flux is sufficiently 
large. Such effect was recently observed in |44 where the 
flux was increased by factor of 100. 



Another deviation from Eq. (4.14) might occur at very 
low flow rate. Here the percolation effect disappears and 
the grains segregate due to size, because the larger grains 
do not find large enough holes in the surface so they roll 
easier than smaller grains. In this case, the rolling phase 
is thin, so that it behaves as a homogeneous phase with 
a constant velocity v fo r al l the grains in the fluid phase. 
In this case, from Eq. (|4.5| ) we obtain 



A = 



Thus when R° 



A*fluid 
Mbulk 

> 



v 



R° 



A - »/(C l7 ) 



(7.2) 



(7.3) 



where d is a small constant of the order of the grain size. 
This lower limit might indicate the crossover from a per- 
colation regime to a non-percolation regime at low flow 
rates. 



In summary, we tested experimentally the main as- 
sumptions of the theory of surface flows of granular ma- 
terials. We measured the profile of the static and rolling 
phases, in order to study the profile of the kink and the 
effects of percolation in the rolling phase and compared 
with analytical predictions. We characterized the dy- 
namical process of stratification by measuring all the rel- 
evant quantities. We measured the velocity of the rolling 
grains, the velocity of the kink, the wavelength of the lay- 
ers, and also the rate of collisions 7, and all four angles 
of repose Q a p characterizing the mixture. The velocity 
of the kink and the wavelength of the layers both vary 
linearly with the grain flux. The velocity profile of the 
grains in the rolling phase is a linear function of the posi- 
tion of the grains along the moving layer, which implies a 
linear relation between the mean velocity and the thick- 
ness of the rolling phase. We find the mean velocity of 
the rolling grains is the same as the velocity of the kink. 

Our results apply to the moderate flow rates used in 
this work of the order of 1 gr/sec. For sufficiently larger 
or smaller flow rates, deviations might appear as dis- 
cussed in the tex t. F or larger fluxes, nonlinear terms 
may modify Eq. ( 4.14 ). For even larger fluxes the kink 
may not support the incoming grains turning stratifica- 
tion into the mixing of grains or to weak segregation. 
For smaller fluxes than the ones used in this study, the 
percolation effect does not take place. However, when 
the size ratio is large enough — p > 1.5 — strong segre- 
gation occurs anyway at the shear surface between the 
fluid and solid phase: the large grains are not trapped 
in the holes of the surface, and they are convected fur- 
ther. Thus stratification is also observed for small fluxes, 
but the size segregation mechanism is different from the 
one studied here. The sharp segregation profiles with a 
separation zone of a few cm observed in the experiment 
shown in Fig. |]b is also a consequence of strong seg- 
regation effects. When p < 1.5, size segregation has a 
weak effect, resulting in a weak continuos segregation of 
the mixture no matter the shapes of the grains. Theo- 
retically the case p < 1.5 is treated in Q, and the case 

p > 1.5 in H-EE}. 

Our results might be also applicable to other systems. 
Size and shape segregation in rotating drums may be 
analized in analogy to the regimes found here. Further 
experimental results may include qualitative studies of 
the other flows regimes mentioned above. It would be 
also appropriate to have an estimation of the angles of 
repose of the grains independently of the theoretical cal- 
culations used here. For instance, by gluing grains of one 
species to an inclined plane and pouring grains of the 
other species is a way to obtain a direct estimation of 
the cross angles of repose. 
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TABLE I. Results of the six experiments. 



Experiment # 


R° (cm) 


A (cm) 


Wf (cm/sec) 


v (cm/sec) 


1 


0.31 


0.63 


6.7 


7.1 


2 


0.38 


0.64 


8.5 


8.7 



TABLE II. Typical values of phenomenological constants. 



7 (1/sec) 


021 


9u 


#22 


012 


Wf (cm/sec) 


v (cm/sec) 


t>/7 (cm) 


23 


8° 


26° 


39° 


57° 


7-18 


7-17 


0.3-0.8 



